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Recent dramatic advances in nanofabrication made it possible to design, arrange 
and assemble nanoparticles with great precision. Optical and, more generally, 
electromagnetic properties of nanostructures have been of great interest in the past 
decade [1-5]. In particular, physical effects due to giant local field enhancement are 
subject of active current research [5-8]. In this respect, nanoparticles of noble metals, 
especially silver, proved to be very useful. The remarkable optical propertied on silver 
nanostructures are explained by the strong, resonant interaction with electromagnetic 
fields in the visible and near-IR spectral range and by very small Ohmic losses. 

The strong enhancement of local fields in small spatial areas is the consequence 
of two factors: the heterogeneity of a nanostructure on a subwavelength scale and the 
resonant character of interaction of the electromagnetic field with the nanostructure. 
Both features are, in principle, present even in the case of a single isolated nanosphere. 
However, the effect becomes much stronger in aggregated nanospheres due to the effect of 
plasmon hybridization [9]. In this case, amplification of local field can become sufficiently 
large to make possible detection of Raman radiation from single molecules, as was 
demonstrated experimentally in [10]. The Raman enhancement factor (|E|/|E |) 4 in 
the center of a junction between two nanospheres (E and Eo - the local and the external 
fields in the junction, respectively) was calculated to be 5.5 • 10 9 for a lnm gap between 
two silver spheres of 60nm radius each at A = 497nm [11]. Even larger enhancement, up 
to 10 13 , was predicted in the so-called nanolens - a linear chain of several nanospheres 
of different size [6] . 

The primary focus of research in single-molecule spectroscopy has been on non- 
coherent optical processes such as Raman scattering. In this paper I consider a coherent 
nonlinear effect, namely, degenerate third order nonlinearity, and demonstrate that, by 
changing the geometry of a nanostructure, it is possible to control not only the amplitude 
of the nonlinear response but also its phase (relative to the phase of the incident field). 
The nonlinear Kerr effect described by the third-order susceptibility x { UJ 'i UJ i UJ i~ UJ ) 
is responsible for corrections to absorption and refraction. The possibility to control 
the phase and tensor structure of x^ 3 ^(o;; to, to, —to) can have numerous applications, for 
example, in quantum nondemolition measurements via the optical Kerr effect [12]. 

Consider a simple physical system shown schematically in Fig. 1. Here a 
nonlinearly-polarizable molecule is placed in the center of symmetry of two spheres. 
The radius of each sphere is denoted by a and the width of the gap (surface-to surface) 
by 2h. Thus, for example, if a = 50nm and 2h = lnm, which is the smallest physical gap 
considered in Ref. [11], we have h/a = 0.01. The system is excited by an external field 
with some frequency to and wavelength in vacuum A = 2-kc/u, linearly polarized along 
the axis of symmetry. The latter coincides with the z-axis. The whole system is assumed 
to be sufficiently small compared to A and we work in the quasistatic approximation. 
Further, we use the Drude formula for the dielectric function of metal, namely, 



e = e - ujI/u(uj + ij) . (1) 
The following parameters are used in the simulations: e = 5, j/tOp = 0.002, and the 
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Figure 1. Schematic illustration of the physical system considered in this paper. 



value of uj p is unspecified. Note that the dielectric function of silver in the anomalous 
dispersion region is well described by choosing uj p « 4.6sec _1 (A p ~ 136nm). Finally, we 
assume that the metal nanoparticles are embedded in a transparent host medium with 
the refraction index = 2 (e^ = 4). 

Because of the axial symmetry, the dipole moment induced in the molecule is 
parallel to the z-axis. The third-order nonlinear correction to the dipole moment 
oscillating at the same temporal frequency as the incident field is given by 

d{ NL \t) = a^E z \E z \ 2 = a<$E \E Q \ 2 exp(-iut) . (2) 

Here E z is the amplitude of local electric field at the location of the molecule and Eq 
is the amplitude of the external field. Note that, even if we assume for simplicity that 
E is purely real, E z can be complex. In general, there can be an arbitrary phase shift 

(3) 

between the local and the external fields. The effective nonlinear polarizability a^ ff is 
related to by 

= Ga^ , (3) 
where the enhancement factor G is given by 
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Since G is, in general, complex, it can influence not only the magnitude but also the 
phase of the effective third-order polarizability. Below, we calculate G numerically 
and show that its phase can be varied in its whole range by changing the inter-sphere 
separation. 

To calculate the local field in the gap, E z , we expand the polarization inside each 
sphere in the quasistatic vector spherical harmonics X^(r) = (/a) -1 / 2 VV^ (r — r,j). 
Here % — 1,2 indexes the nanospheres, are the radius- vectors of the spheres' centers, 
Wm( r ) = ( r / a )'^m(f), an d Yi m (r) are spherical functions of the polar angles of the unit 
vector r. Polarization inside each sphere can be written as 

P(r) = CflnvXSLOO , if |r - r 4 | < a , (5) 
i in 

where the unknown coefficients Cu m must be found from the standard boundary 
conditions applied at the surface of each sphere, or alternatively, from the integral 
equation formalism as described in [13]. From general considerations, it is clear that 
Cum obey a system of linear equations which, in the quasistatic limit, was obtained 
in [14] and simplified in [15]. In general, this set of equations has the form 

(1/ X -W)\C} = \E) , (6) 

where x — (3/47r)[(e — £h)/( £ + % £ h)] is the coupling constant, W the electromagnetic 
interaction matrix, and \E) is the appropriate right-hand side defined by the external 
field. In the case of axial symmetry, only modes with m = are excited, so that 
Cam = Cu 5 m0 . The matrix elements of W needed to find the solution are 



(M\W\i'l>0) = + (1 - fe)(-lf [sgn(z, - z if )\ l+v 



(21 + 1)(2/' + 1) + h/a) l+l ' +1 

and the components of the right-hand side vector in (6) are given by 

(M\E) = E ^/4:7ia 3 /3 . (8) 

In general, once the coefficients Cu m are found, the scattered field at an arbitrary 
point r in the host medium can be found from 

E s (r) = - J2 C *m S^y X !?i( r ) > (9) 

where X^(r) = [(/ + l)a] _1 / 2 V^^(r — r^) are the quasistatic vector spherical harmonics 
of the second kind and 4>hl{ r ) = ( a / r )' +1 ^»n( r )- For the particular problem considered 
in this paper, summation (9) contains only terms with m = 0. Further simplification is 
obtained if the electric field is evaluated on the axis of symmetry, in which case 
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where 



/„(*) = 



oo 



(-1) 



(2/0|P n ) 




(11) 



z=l 



2Z + 1 \[l + (/t/a)(l+i)] 



Z+2 



[l + (Va)(l-x)] 



and |P n ) are the eigenvectors of W with corresponding eigenvalues w„. Here E sz is the 
z-component of the scattered field on the axis of symmetry (the x- and y-components 
are zero). The total local field E z is a superposition of the incident and scattered fields: 



Note that we have used the spectral approach to solving (6). In other words, instead of 
directly inverting 1/x — W, we seek eigenvectors and eigenvalues of W and then obtain 
the solution in terms of these quantities for an arbitrary coupling constant x- 

The matrix W is of infinite size and in practical calculations must be truncated. 
The truncation order Z max required to obtain an accurate solution depends on the inter- 
sphere separation. Although, for any separation, there exist an infinite number of modes, 
most of them are antisymmetric [16,17] or, equivalently, dark [18]. A dark mode is 
not coupled to the homogeneous external field because the scalar product (E\P n ) is 
either exactly zero or very small. Correspondingly, the input of a dark mode to sum 
(10) is negligible. Modes which are not dark are referred to as luminous [18]. For 
a finite value of h, there is a finite number of luminous modes and the spectrum of 
eigenvalues w n which correspond to these modes is discrete. However, as h decreases, 
the number of luminous modes grows and the intervals between consecutive values of 
corresponding eigenvalues w n approach zero. When the two spheres touch, the spectrum 
becomes continuous and, strictly speaking, can not be adequately described at any finite 
truncation order. However, for practical purposes, the matrix still can be truncated, 
as long as the resultant discrete density of states approximates the true continuous 
function with sufficient precision. The latter condition depends on the relaxation in the 
system and is very difficult to satisfy for silver in the near-IR spectral region due to the 
very small non-radiative relaxation. In the simulations presented below, the minimum 
ratio h/a is equal to 0.01. In this case, all luminous modes are obtained with very 
high precision at relatively modest truncation orders. The results reported below were 
obtained at / max = 800 and convergence with machine accuracy was verified by doubling 
this number. 

Now we turn to the numerical results. First, in Fig. 2, we plot the spectral 
dependence of the factor \E z /Eq\ 2 in the center of the inter-sphere gap for different 
relative separations h/a. It can be seen that a "resonance band" exists in the spectral 
region whose bounds depend on the ratio h/a. For h = 0.01a, resonance interaction 
takes place for 0.15 < uj/uj p < 0.33. For h = 0.32a, the resonance band is smaller, 
0.25 < uj/uj p < 0.32. We will be interested in the frequencies which lie in the resonance 
band for the smallest value of h considered, namely, h = 0.01a. 



E z — E + E, 



(12) 
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Figure 2. The ratio \E z /Eq\ 2 in the center of the inter-sphere gap as a function of 
Lu/ujp for different relative inter-sphere separations. 

In Figs. 3-5, we show the parametric plots of the complex enhancement factor G for 
the following values of the ratio u/oj p . 0.20, 0.25 and 0.32. In the case u/u p = 0.20, the 
most dramatic change of G happens when h changes from 0.36a to 0.41a. The phase 
of G changes in this interval of h from pa it/ '4 to pa 37r/4. Overall, the phase of the 
enhancement factor can be tuned from pa to pa ir by tuning h in the whole considered 
interval. 

Much more control over the phase of G can be attained for u/u p = 0.25, as shown 
in Fig. 4. It is interesting to note that the parametric curve shown in these figure is 
approximately self-similar, consisting of several almost closed loops which can be seen at 
different scales. The phase of G changes in the whole interval from to 2%. Qualitatively 
similar curve was also obtained for uo/uo p = 0.30 (data not shown). 

Perhaps, the most interesting curve is obtained at uj/uj p = 0.32 (Fig. 5), although 
the magnitude of G is not as large for this value of u/u p as in Figs. 3-4. The parametric 
plot of G is in this case a spiral. The phase of G changes monotonously from pa 
to pa 57r. Thus, the curve makes more than two full revolutions around the origin in 
the complex plane. By varying the parameters in the Drude formula and the refractive 
index of the host medium, it was found that the spiral shape of the curve is typical when 
uj/uj p is close to the right bound of the resonance interaction band (data not shown). 

Thus, we have shown that by changing the inter-sphere separation h from 0.01a to 
0.32a it is possible to change the phase of the enhancement factor G, and, consecutively, 
that of the effective nonlinear polarizability a e{ { in its whole range. Some limitations 
of the model used in this paper must be mentioned. First, we did not account for 
direct electromagnetic coupling between the nonlinearly polarizable molecule and the 
nanospheres, nor did we take into account the nonlinearity of the metal itself. The latter 
effect can be significant. Further, did not account for the fact that the electric field in 
the gap is not constant but can change on the scales comparable to molecular. Although 
all these factors are important if one seeks to calculate the nonlinear response of the 
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Figure 3. Parametric plot of the complex enhancement factor G as a function oih/a 
for uj/ui p = 0.2. Graphs (a)-(c) show the same curve on different scales. 
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Figure 4. Same as in Fig. 3 but for to = 0.25u> p . 
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Figure 5. Same as in Fig. 3 but for u = 0.32o;p. 

system with precision, inclusion of all these complications would make the theoretical 
description intangible. On the other hand, the physical effect described in this paper 
does not originate due to any of the approximations listed above. Instead, it is explained 
by the resonant nature of the interaction between the electromagnetic field and the 
nanosystem. When the spacing between the nanospheres is tuned, different resonance 
modes are excited in the bisphere aggregate. This results the change of the relative phase 
between the local field in the gap E z and the external field Eq and the characteristic 
dependence of the enhancement factor G on h which is illustrated in Figs. 3-5. Since 
the resonance nature of interaction is not altered by the effects mentioned above, it 
is reasonable to expect that the fine tuning of the third-order nonlinear response is 
achievable in nanosystems specifically engineered for that purpose. 
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